Rapid granular flows on a rough incline: phase diagram, gas transition, and effects of air drag

We report experiments on the overall phase diagram of granular flows on an incline with emphasis on high inclination angles where the mean layer velocity approaches the terminal velocity of a single particle free falling in air. The granular flow was characterized by measurements of the surface velocity, the average layer height, and the mean density of the layer as functions of the hopper opening, the plane inclination angle and the downstream distance x of the flow. At high inclination angles the flow does not reach an x-invariant steady state over the length of the inclined plane. For low volume flow rates, a transition was detected between dense and very dilute (gas) flow regimes. We show using a vacuum flow channel that air did not qualitatively change the phase diagram and did not quantitatively modify mean flow velocities of the granular layer except for small changes in the very dilute gas-like phase.Comment: 10 pages, 16 figures, accepted to Phys. Rev.

Slippery sliding on icy Iapetus

Enigmatically, some landslides flow farther than normal frictional resistance allows. Cassini images of Saturn’s icy moon Iapetus reveal a multitude of long-runout landslides that may have been enabled by flash heating along the sliding surface

From hot rocks to glowing avalanches: Numerical modelling of gravity-induced pyroclastic density currents and hazard maps at the Stromboli volcano (Italy)

Gravity-induced pyroclastic density currents (PDCs) can be produced by the collapse of volcanic crater rims or due to the gravitational instability of materials deposited in proximal areas during explosive activity. These types of PDCs, which are also known as “glowing avalanches”, have been directly observed, and their deposits have been widely identified on the flanks of several volcanoes that are fed by mafic to intermediatemagmas. In this research, the suitability of landslide numerical models for simulating gravity-induced PDCs to provide hazard assessmentswas tested. This work also presents the results of a back-analysis of three events that occurred in 1906, 1930 and 1944 at the Stromboli volcano by applying a depth-averaged 3Dnumerical code named DAN-3D. The model assumes a frictional internal rheology and a variable basal rheology (i.e., frictional, Voellmy and plastic). The numerical modelling was able to reproduce the gravity-induced PDCs' extension and deposit thicknesses to an order ofmagnitude of that reported in the literature. The best resultswhen comparedwith field datawere obtained using a Voellmymodelwith a frictional coefficient of f=0.19 and a turbulence parameter ξ=1000 m s−1. The results highlight the suitability of this numerical code,which is generally used for landslides, to reproduce the destructive potential of these events in volcanic environments and to obtain information on hazards connected with explosive-related, mass-wasting phenomena in Stromboli Island and at volcanic systems characterized by similar phenomena.Published93-1065V. Dinamica dei processi eruttivi e post-eruttiviJCR Journa

New, simplified and improved interpretation of the Vaiont landslide mechanics

Both the occurrence and behaviour of the Vaiont landslide have not been satisfactorily explained previously because of difficulties arising from the assumption that the failure surface was ‘chair’ shaped. It is now known that there was no ‘chair’, which means that the 1963 landslide could not have been a reactivated ancient landslide because the residual strength of the clay interbeds would have been insufficient for stability prior to 1963. Furthermore, the moderately translational geometry reduces the influence of reservoir-induced groundwater and hence of submergence. Standard stability analyses now show that prior to 1960, the average shear strength must have significantly exceeded the peak shear strength of the clay interbeds known to have formed the majority of the failure surface. Three-dimensional stability analyses confirm these results and show that at the time of the first significant movements in 1960, the rising reservoir level had a negligible effect on the Factor of Safety. According to these results, the Vaiont landslide was most likely initiated by pore water pressures associated with transient rainfall-induced ‘perched’ groundwater above the clay layers, in combination with a smaller than hitherto assumed effect of reservoir impounding, then developed by brittle crack propagation within the clay beds, thus displaying progressive failure. Further, very heavy rainfall accelerated the process, possibly due to reservoir-induced groundwater impeding drainage of the rainwater, until the limestone beds at the northeast margin failed. With the shear strength suddenly reduced to residual throughout, the entire mass was released and was able to accelerate as observed

Coriolis-induced instabilities in centrifuge modeling of granular flow

Granular flows are typically studied in laboratory flumes based on common similarity scaling, which create stress fields that only roughly approximate field conditions. The geotechnical centrifuge produces stress conditions that are closer to those observed in the field, but steady conditions can be hardly achieved. Moreover, secondary effects induced by the apparent Coriolis acceleration, which can either dilate or compress the flow, often obscure scaling. This work aims at studying a set of numerical experiments where the effects of the Coriolis acceleration are measured and analyzed. Three flow states are observed: dense, dilute, and unstable. It is found that flows generated under the influence of dilative Coriolis accelerations are likely to become unstable. Nevertheless, a steady dense flow can still be obtained if a large centrifuge is used. A parametric group is proposed to predict the insurgence of instabilities; this parameter can guide experimental designs and could help to avoid damage to the experimental apparatus and model instrumentation

Identification of debris flow initiation zones using topographic model and airborne laser scanning data

Empirical multivariate predictive models represent an important tool to estimate debris flow initiation areas. Most of the approaches used in modelling debris flows propagation and deposit phases required identifying release (starting point) area or source area. Initiation areas offer a good overview to point out where field investigation should be conducted to establish a detailed hazard map. These zones, usually, are arbitrarily chosen which affect the model outputs; hence, there is a need to have accurate and automated means of identifying the release area. In addition to this, the resolution of the terrain dataset also affects the results of the detection of source areas. In this study, airborne laser scanning (ALS) data was used because of its robustness in providing detailed terrain attributes at high resolution. Primary and secondary conditioning parameters were derived from digital elevation model (DEM) as input into the modelling process. Three models were executed at different spatial resolution scales: 5, 10 and 15 m, respectively. MARSpline multivariate data mining predictive approach was implemented using morphometric indices and topographical derived parameter as independent variables. A statistics validation was calculated to estimate the optimal pixel size, 1200 randomly sample data were generated from existing inventory data. Debris flows and no-debris flows were categorized, and the transform to continuous integer (1 and 0), respectively. To achieve this, the data set was divided into two, 70% (840) for the training dataset and 30% (360) for validation. The best model was selected based on the model performance using the generalized cross validation (GCV) and the receiver operating characteristic (ROC) curve/area under curve (AUC) values. Conditioning parameters were numerically optimized to identify the arbitrarily maximum model basis function for eleven variables, using MARSplines analysis (algorithm). The three most influencing topographic parameters identified are topographic roughness index (TRI), slope angle, and specific catchment area (SCA) with the percentage values of participation in the model of 100, 93, and 86%, respectively. The chosen function appeared to describe the analysed correlation sufficiently well. Consequently, three stages of optimization were made to determine the optimized source areas is possible with 10 m pixel size, 200 maximum basis functions and 3 maximum interactions, resulting into 82% ROC train and 80% test, GCV 0.189 and 85% correlation coefficient. The result will be of great contribution to the advancement of a broad understanding of the dynamics of debris flows hazard and mitigations at regional level which; that is resourceful for comprehensive slope management for safe urban planning decision-making process and debris flow disaster management